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Abstract 

We investigate non-stationary orthogonal wavelets based on a non-stationary 
interpolatory subdivision scheme reproducing a given set of exponentials. The con- 
struction is analogous to the construction of Daubechies wavelets using the subdi- 
vision scheme of Deslauriers-Dubuc. The main result is the smoothness of these 
Daubechies type wavelets. 
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1 Introduction 

The main purpose of the present paper is to study a non- stationary interpolatory subdi- 
vision scheme reproducing exponentials and the related Daubechies type wavelets. 

In [9] G. Deslauriers and S. Dubuc investigated a subdivision scheme based on polyno- 
mial interpolation of odd degree 2n — 1. They proved the existence of the basic limit func- 
tion $^271 ^Qj. ^j^g subdivision scheme and determined its order of regularity. Daubechies 
wavelets are closely related to this construction since the autocorrelation function of the 
Daubechies scaling function is equal to the basic limit function $^2n (^ggg jgj^ |33j^ qj. ggg_ 
tion 5 below). Moreover, the regularity of the Daubechies wavelets can be derived from 
the regularity of the basic limit function $^2™ 

The motivation for the present work originated from the attempt to define a new con- 
cept of multivariate subdivision scheme in the spirit of Deslauriers and Dubuc, which is 
based respectively on multivariate interpolation. According to the Polyharmonic Paradigm 
introduced in [19] polyharmonic interpolation on parallel hyperplanes (or on concentric 
spheres) provides a proper generalization of the one- dimensional polynomial interpolation 
and, respectively, generates a natural multivariate subdivision scheme. We shall present 
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in a forthcoming paper [T3] the polyharmonic subdivision scheme as an immediate general- 
ization of the one-dimensional subdivision scheme of Dubuc-Deslauriers. This subdivision 
scheme is stationary and reproduces polyharmonic functions of fixed order. It is reduced 
to an infinite family of one-dimensional inherently non- stationary subdivision schemes 
reproducing exponential polynomials of special type. Thus the results about subdivision 
schemes for exponentials are the building bricks of the polyharmonic subdivision schemes 
discussed in detail in [13], [20], [21] • 

On the other hand, the study of non-stationary subdivision schemes reproducing 
general exponential polynomials were initiated in 2003 in [H], where the main results of 
Deslauriers and Dubuc were generalized, in particular, existence of basic limit function 
and its regularity were proved. For a given set of numbers {Aj}"^q , called sometimes 
frequencies, such a scheme reproduces the space V = span |e'^-'*}"_g, and is character- 
ized by a family of symbols {a^''^{z)}kez, where a^''^{z) defines the refinement rule at level 
k. Ch. Micchelli proved in [30] that the symbols of these schemes are non-negative on 
the unit circle, whenever the set of frequencies {Aj}"^Q are real and symmetric. Based 
on this property of the symbols, he applied in [30j the construction of Daubechies to 
the fixed symbol a^^^{z), corresponding to a set of exponentials defined by the frequen- 
cies {Ai, A2, A„, — Ai, — A2, — A„, 0} . However, the so constructed wavelets reproduce 
only the function e° = 1. In contrast to that, in the present paper we elaborate on a 
genuine non- stationary scheme: we apply the same construction to the family of sym- 
bols {'^''^K^;)}^^^ of a non-stationaly scheme reproducing V, and show that the resulting 
non-stationary wavelets generate a Multiresolution Analysis (MRA) which contains V . 
Let us mention that in [36] a similar construction is proposed without the validation of 
the positivity of the symbols for all k, which is necessary for the construction; see more 
discussion on the above references at the end of the paper. 

The main result of the present paper states that the order of regularity of the new 
Daubechies type wavelets reproducing e^^^ for j = 1, ...,n, is at least as large as in the 
case of classical Daubechies wavelets reproducing polynomials of degree <n — l provided 
that the filters of the Daubechies scheme are chosen at each level in an appropriate way. 
The proof of this result depends on the concept of asymptotically equivalent subdivision 
schemes developed in [16]. These regularity results are important for the regularity of 
the multivariate polyharmonic subdivision schemes and the corresponding multivariate 
wavelets considered in [T5] . 

What concerns the computational aspects of the present results, let us mention that 
in the case of a general set of frequencies {A-, }"^-,^ it is difficult to find concise expressions 
for the symbols of subdivision and respectively for the filters of the wavelets, in particular 
for the analogs of the important classical polynomials Qn (appearing in equality (l25l) 
below), and the polynomials U'^^z) (in equality (1321) below). However we have to note 
the remarkable fact that in the special cases of interest for the multivariate polyharmonic 
subdivision scheme the corresponding polynomials can be explicitly constructed, see [13j, 
[20] . [21], [22], [23]. Also, in [20j, [21], [^, [23] the polyharmonic subdivision wavelets 
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of Daubechies type were successfully applied to problems of Image Processing, where 
the experiments were carried out on benchmark images as well as on some astronomical 
images. 

Finally, let us mention that after jSjj non-stationary Wavelet Analysis reproducing 
exponentials have appeared in a natural way in the monograph [19] in the context of 
multivariate polyspline wavelets. 

The paper is organized as follows: in Section 2 we shall review basic notions for 
non-stationary subdivision schemes. Section 3 is devoted to subdivision schemes for ex- 
ponential polynomials and we shall give an improved estimate of the order of regularity 
of the basic limit function. In Section 4 we shall briefly review non-stationary multires- 
olutional analysis. In Section 5 we shall give present the construction of the Daubechies 
scaling function via subdivision schemes and estimate the order of regularity of the new 
Daubechies type wavelets. 

Finally, we introduce some notations: No denotes the set of all natural numbers includ- 
ing zero, Z denotes the set of all integers and we define the grid 2~^Z = : j G Z} . 
By (R) we denote the set of all functions / : M — > C which are £-times continuosly 
different iable where £ G Nq. The Fourier transform of an integrable function / : M — ?■ C is 
defined by 

oo 

f{uj) = I f{x) e-'^'^dx. 



2 Basics in nonstationary subdivision schemes 

In this section we shall briefly recall notations and definitions used in non- stationary 
subdivision schemes for functions on the real line. The formal definition of a subdivision 
scheme is the following: 

Definition 1 A (non-stationary) subdivision scheme Sq is given by a family of se- 
quences (a}''^] of complex numbers indexed by k G Nq, called the masks or rules at 

level k, such that a^^^ ^ only for finitely many j G Z. Given a sequence of numbers 
/° (j) ) J ^ ^) one defines inductively a sequence of functions f^^^ : 2~^^~^^^'Z C by the 
rule 

fr (^) = (i) /"'-^^ 2- 

If for each /c G No the masks € Z, are identical the scheme is said to be stationary. 
The subdivision scheme is called interpolatory if for all k>0 and j G Z holds 

f2^' = ft (2) 
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An important tool in subdivision schemes is the symbol a'*^^ of the rule, or mask 
G Z, defined by 

\k] 

which is Laurent polynomial since we assume that the sequence uj ,j G Z, has finite 
support. We shall identify the subdivision scheme 5*0 by its masks af\j G Z, G No or 
its symbols a^''\ k G Nq. It is clear that a subdivision scheme is stationary if and only if 
a^'^^lz) = a{z), for all k G Nq. Moreover the scheme is interpolatory if and only if afj = 6oj, 
for all j G Z. In terms of the symbol a'*^' it is easy to prove that this is equivalent to the 
identity 

al'=l(z) + aW(-z) = 2, ^ G C \ {0}. 

Definition 2 Let i G Nq. A subdivision scheme So is called -convergent if for any 
bounded initial sequence {fj}jez there exists F E (M) such that 

lim sup { If (j2-'') - /m = 0. (3) 

The function F is called limit function of the subdivision scheme for the initial 
sequence {fj}j(zz- The limit function for the initial data function fj = S^j (here Sqj is the 
Kronecker symbol) is called basic limit function of the scheme and it is denoted by $o- 

A subdivision scheme is called convergent if it is just C^'-convergent. A central 
problem in the theory of subdivision schemes is to estimate the order of regularity of a 
limit function whenever it exists. Let us remark that for an interpolatory scheme the 
limit function F in Definition [2] is easily computed for dyadic numbers t = ^ with k > 
and j E 1j through the values f^ and convergence of the scheme asks whether it has a 
continuous extension over the whole R. 

Remark 3 It is easy to see that for a convergent subdivision scheme Sq with symbols 
_,^af for all k G No, the support of the basic limit function $o is 
contained in [—N, N] . 

In the non-stationary case the following concept is of importance: 

Definition 4 Let Sq be a subdivision scheme given by the masks (af^) . For any 

natural number m G Nq define a new subdivision scheme Sm by means of the masks of 
level k : 

flW'™ (z) := a['=+'"l (z) for all k G Nq. 
The following result is proved in [15j : 
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Theorem 5 // 5*0 is a convergent subdivision scheme then Sm is convergent for any m G 
No. 

The basic limit function of this scheme is denoted by 

The foUowing resuh is well-known, see e.g. the proof of Theorem 2.1 in [5], for further 
results see also |311. 



Proposition 6 Let Sq be a subdivision scheme with symbols a^''^{z) for k eNq such that 
(i) a^'^' = for all \j\ > N and all G Nq for some fixed integer N > 0. 
(a) there is a constant M > such that lal'^l < M for all /c G Nq and 



E 

A;=0 



-aW (1) - 1 
2 ^ ' 



< oo. (4) 



Then the infinite product 



2" 



k=l 

converges uniformly on compact subsets o/M. 

Note that by Remark 3 the basic limit function $o of ^ convergent subdivision scheme 
satisfying (i) in Proposition |6] has compact support. Hence the continuous function $o 
is integrable and square integrable. It follows that the Fourier transform $o of $o is 
well-defined, continuous and square integrable. 

Proposition 7 Let Sq be a convergent subdivision scheme with symbols a^'^^z) for k E Nq 
satisfying (i) and (ii) in Proposition\^ Then 

oo ^ 



k=l 



\k] 

Definition 8 A subdivision schemes Sq with masks aj ,j E Z, k E No reproduces a 
continuous function / : M — t- C at level k eNq if for all j E 



2fc+i 

We say that Sq reproduces f stepwise if it reproduces f at each level k E Nq. 

In [151 p. 31, just called reproducing] this is called just reproducing, and in [121 here 
"stepwise"] is defined the stepwise reproduction. 
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Proposition 9 Let Sq be a subdivision scheme with symbols a^'^^z) for k G Mq, let A G C, 
and g (x) be a polynomial. Define z^'^'^ := exp (_2-('=+i)A), and for 6 = 0,1 define 



for 6 = 0,1. 

Then for fixed /c G Nq ^/ie following statements are equivalent : 

a) g (x) e^^ is reproduced at level k. 

b) Ffl [g, m)=g (|g±^) for 5 = 0, 1 and for all m G Z. 

c) fP {g, x) = g (f|#) for 5 = 0,1 and for all x G M. 

Proof. Using ([5]) for j = 2m + 5 with m G Z and 5 = 0, 1, it is easy to see that the 
function g {x) e^^ is reproduced stepwise if and only if for 5 = 0, 1 and for all m & Z 



I 2'=+! I ~ ^ ^2m+5-2l9 1 2^ / ^ ^ 



2m+S-2l 



Using the variable transformation / = / + m it is easy to see that the right hand side is 
equal to Fs {g,rn) , hence the equivalence of a) and b) is proven. The equivalence of b) 
and c) follows from the fact that g (x) and Fs {g, x) are polynomials. ■ 

Theorem [10] below extends Theorem 2.3 in [13] which was formulated only for in- 
terpolatory schemes. (We will apply this later to Daubechies schemes which are not 
interpolatory) . 

Theorem 10 Let Sq be a subdivision scheme with masks cif\j G Z, A; G Nq. Then for 
r = 0, — 1 the functions fr (x) = x'^e^^ are reproduced stepwise by Sq if and only if 
for all A; G No holds 

(- exp (-2-('=+i)A)) = and (exp {-2-^''+^^)) = 2 (6) 

and 

|^aW(±exp(-2-(^-+^)A))=0, r = 1, ...,/.- 1. (7) 

Proof. Put := exp (-2"('=+i)A) . Since igs,0) = Eiez^Lds (^) {z^'^Y'"' , 
we first observe the following identities for the constant function g = 1 : 

aW(^W) =fJ'=] (1,0) + Ff 1(1,0), (8) 
aW (-^W) = Flt^ (1, 0) - Ff 1 (1, 0) . (9) 

Suppose 5*0 reproduces the function f (x) = g (x) e^^ where g (x) is a polynomial of degree 
< /i — 1. We obtain from condition b) in Proposition [9] for g = 1 that 1 = Ff^^ (1,0) for 
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(5 = 0, 1, and we conclude from and © that (^W) = 2 and (-zl*^') =0. Next 
we take 

gl (x) = {6- {6 - 2''+^x - 1) ... {6 - - {s - 1)) (10) 

for s = 1, /i — 1. Then 

fil = ('^ - 2/) ('^ - 2/ - 1) ... (5 - 2/ - (s - 1)) . (11) 



2k ^ 

Since Ff (^75,0) = Ez6z4l2^^?<5 (2^) {z^'^Y''' it follows from (HI]) that 

^aW(.W)=i^f (^?o,0) + Ffl (,^,0), (12) 



-a 



W =FW(^^,0)-Ff](^,^0). (13) 



Since (7I (x) e^^ is reproduced, we obtain that f]'^^ ((/I, 0) = gs (2^) = for 5 = 0, 1. Thus 
([7]) holds and the necessity part is proved. 

The converse is proved by induction over the dimension s of the linear space Eg of 
all polynomials g of degree < s. By Proposition [9] we have to prove that for each g ^ Eg 
holds 

9 ( = i9,m) for 5 = 0,1 and for all m G Z. 



2k+i 

For s = this means we have to prove that Ff'^ (1,?^) = 1 for the constant function 
g = 1, for 5 = 0, 1, and for all m G Z. Since f]'^' {g, x) is the constant polynomial we have 
only to show that f]'^' (1,0) = 1. By assumption ([7]) and equations ([8]) and (Q we infer 
2 = aW (zW) = Fj'^l (1, 0) + Ffl (1, 0) and = = Fj'^' (1, 0) - Ff ^ (1, 0) , from 

which we infer the desired statement f]'^^ (1, 0) = 1 for 6 = 0,1. 

Suppose that the statement is proven for Eg-i- We consider the polynomial g^ defined 
in (fTOl) for 5 = 0, 1. Note that equation (fT2|) and (|T3|) . and the assumption ([7j) imply that 

Fl'\9hO)=0 = gl(^-^^ (14) 

for all s = 1, ...,fi- 1. Put now g' := ^l^. Then by it follows F^''' (^^ 0) = g' (0) . 
Note that 

gl (x) = g^, (x) + G^-' (x) 
for some polynomial G^"^ of degree < s — 1. Then 



Fl'^ {g'o, 0) = Ff] (^^ 0) - Ff 0) = -G-^ 



2fc+i 
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since Ff ' (gl 0) = by ([HD, and ' {G'-\ 0) = G''^ {^) by the induction hypothesis. 
Hence 
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2fe+i 



If we prove that fJ^^ {g% m) = g'^ (^ftt) for allm G Z and 5 = 0, 1, then we can conclude 
that g"^ is reproduced by the scheme and clearly this implies that Eg is reproduced. Note 
that for each m G Z there exists a polynomial of degree < s such that 

9' + ^) = 9' (x) + /i^ (x) . 

Then + f ) = g' {^) + (^) and 

By the induction hypothesis we know that f]'^^ {h^, 0) = hf^ iw^) ■ follows from f|T^ 
that 

The proof is complete. ■ 

It is easy to see that a function / : M — t- C is reproduced stepwise by a subdivision 
scheme if and only if for the data function /^^l (j) := / (j) one has 



J ' ^2'' / 

From this it is easy to see that a convergent and stepwise reproducing subdivision 
scheme is reproducing in the following sense: 

Definition 11 A convergent subdivision scheme Sq with masks G G Nq, is 

reproducing a continuous function f : ^ C if the limit function for the data function 
fU), 3 ^ ^, is equal to f (x) . 

Remark 12 In the definition of a convergent scheme some authors require bounded data. 

An important concept for the investigation of non-stationary subdivision schemes is 
the following notion introduced in 



Definition 13 Two subdivision schemes Sa and Sb with masks afKj G Z, A; G No, and 

\k] 

by,jGlj,kG No,resp. are called asymptotically equivalent if 

oo 

^^|af-(,f|<«3. (15) 

fc=0 jGZ 



We say that Sa and Sb are exponentially asymptotically equivalent if there exists a 
constant C > such that 

max 

for all G No- 

We also say that the masks a^'^l, k G No, and b^''\ k E Nq, are (exponentially resp.) 
asymptotically equivalent if ffT^ ( ffTB]) resp.) holds. 

Suppose that the masks 0'f\j G Z, A; G Nq, and b^j^\j G Z, A; G Nq, have support in 
the set {-A^, A^} , i.e. = for \j\ > N. Then it is easy to see that a^''\ k G Nq, 
and fot'^l, G No, are exponentially asymptotically equivalent if and only if for any R > 1 
there exists D > such that 

I^W (z)-6W {z) \ <D-2-'' 

for all A; G No and for all z G C with 1/R< \z\ < R. 

The following result provides a sufficient method for constructing asymptotically equiv- 
alent subdivision schemes. 

Theorem 14 Let m G No and assume that p^^' (z) and p (z) are polynomials of degree m 
for each A; G No defined as 

m m 

pl'^l (z) = n (-^ ~ p{z)=cY[{z- aj) 

i=i i=i 

for some complex numbers c^''\c, and a^^^ and aj, for j = l...,m, and k G No. Suppose 
that there exists a constant D.^ > such that for all A; G No and j = 1, m 

< D^2-^ and jc'^'l - c| < 0^2'''. (17) 

Then p^'^^ (z) ,k G No, and p{z) , A; G No, are exponentially asymptotically equivalent. 

Proof. We claim by induction over m G Nq that for each R> there exists a constant 
Cm {R) such that for all A; G No and \z\ < R 

|/'U^)-p(^)| <C„(r)2-^ (18) 

For m = the statement follows directly from f|T7|) . Suppose that the statement is true 
for m — 1. Write 

p[k] _ ^[k] _ ^[fe] j p{z)=c{z- am) Pm-1 {Z) , 



[k] _ Ak] 
j j 



<C -2- 



(16) 



a 



[k] 



— a. 
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where p^-i (^) Pm~i (z) are polynomials of degree < m — 1 and leading coeffi- 
cient 1. By the induction hypothesis, for each i? > a constant Cm {R) such that 

< Cm (r) 2-^= for all A; G No and \z\ < R. Note that 



{z) - CPm-l {z) 



{Z) -P{Z)= [Z - al^l) {z) -c{z- am) Prn-l (z) 

= {Z- am) [cWpS_i {Z) - CPm-l] + - «W) p^^_, (z) . 

Using the triangle inequality, the induction hypothesis, and ffTTI) . one obtains ffTSl) . ■ 

Theorem 15 Assume that p^''^ (z) and p (z) are polynomials of degree < m for each 
/c G No and assume for each R > there exists a constant Cm (R) such that for all A; G No 
and \z\ < R 

\p^'Hz)-piz)\<CmiR)2~\ (19) 

If a is a simple zero of p (z) then there exists kQ G No and a constant p > such that for 
each natural number k > kQ there exists a zero at'^' of p^''^ (z) with 

\a^''^ -a\< p2-^ for all k > ko- 

Proof. Let p > and define 'jk (t) = a + 2~'^pe** for t G [0, 27r] . Write p (z) = 
Ylf=iPj {z — aY . Since p (z) has a simple zero in a we know that po = and pi ^ 0. We 
obtain the estimate 



b(7fc(t))|>p2-M bl|-5^b,|- 



P^ 



2fc(j--i) 

Hence for given p > there exists a natural number kp such that for all k > kp 

|p(7.W)l>p2-^'^ 

and for all t G [0, 27t] . Now take i? > is large enough that i? > 2 |a| . Take p such that 
p^-^ > Cm (R) where Cm (R) is the constant in (fT9l) . Take fci > kp large enough that 
hk {t)\<R for all A; > A;i and t G [0, 2n] . Then 

b''' ilk (t)) - P {lu {t))\ < Cm (R) 2-' < p^2-'= < \p (7fc m 

for all k > ki and t E [0, 2%] . By Rouche's theorem the number of zeros of p'^^ (z) in the 
ball 1^ — a| < 2~'^p is equal to the number of zeros of p (z) in that ball. Since p (a) = it 
follows that for each k > ki there exists a zero a^'^^ of p^'^^ (z) in the ball |2; — a| < 2~^p, 
hence, la'^'^^ — al < 2"^p. ■ 
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3 Subdivision schemes based on exponential interpolation and 
regularity of the basic limit function 

Let us first recall some basic facts about the classical 2n-point Deslauriers-Dubuc subdi- 
vision scheme. 

It is defined via interpolation of polynomials of degree 2n — 1, see [9]. We shall denote 
its symbol by D2n {z) which is given by 

D2n{z)= Yl PO^'- (20) 

lil<2n-l 

According to Theorem 6.1 in ^ the symbol D2n{z) has the (reproduction) property if 
and only if 

—D2n{-1) = 0, for 3 = 0, 2n - 1. (21) 

Further the scheme of Deslauriers and Dubuc is interpolatory, i.e., p2j = ^oj, for all j G Z, 
or equivalent ly 

D2n{z) + D2n{-z) = 2, for all zeC\ {0}. (22) 

Together, conditions (I2T1) and ( l22l) constitute a linear system which uniquely determines 
the symbol D2n{z) of the Deslauriers and Dubuc scheme and it can be written in the form 

D2n{z)=(^-^X\D,A^)- (23) 



Let us mention that condition ( ]2T]) means that polynomials of degree < 2n — 1 are 
reproduced by the subdivision scheme. The Laurent polynomial 6D2n (z) can be computed 
explicitly and the following identity (see e.g. [T3]). 

D2n (z) = 2^ ''^l g„_i ((^ (z)) (24) 

holds for all -2 7^ 0, where (p (z) = ^ — \{z + l/z) and Qn~i {x) is the polynomial of degree 
n — 1 given by 

The next result shows that Qn-i {z) and fo/jj^ {z) have only simple zeros. 
Proposition 16 The polynomial Qn-i (x) in l[25\) satisfies the identity 



nQn-i (x) + (x - 1) (x) 



(2n- 1)! 



■n-l 



(n - 1)\ (n - r 



tX 
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Proof. Clearly Q[^_^ (x) = ■ ^hen nQn-i + - 1) Q^.^ (x) is equal 



to 

n— 1 / , . n—l , , . n— 2 



from which the statement follows. ■ 

Now we turn to subdivision schemes for exponential polynomials. Let L be the linear 
differential operator given by 

AoV... f:^-A„V (26) 



dx J \dx 

Complex-valued solutions / of the equation Lf = are called L-polynomials or exponen- 
tial polynomials or just exponentials. We shall denote the set of all solutions of Lf = 
by E(Ao...,A„). 

Throughout this article we shall assume that Aq, A„ are real numbers. This assump- 
tion implies that E {Xq..., Xn) is an extended Chebyshev space, in particular for any 
pairwise distinct points to, ...,tn and data values yo, ...,yn, there exists a unique element 
p G E {Xq..., Xn) with p(tj) = fjitj) for j = 0, ...,n, i.e. p is interpolating exponential 
polynomial, cf p6j. 

Given real numbers Aq, X2n-i one can define the subdivision scheme based on inter- 
polation in E (Aq..., A2n-i): the new value /'^^^ (j72'^+^) is computed by constructing the 
unique function pj E E (Aq..., A2ri-i) interpolating the previous data /'^ ((j + /) /2'^) for 
/ = -n + 1, ...,n, and putting (j/2'=+^) = p^ (j/2'=+i) (see [30] for details). Then the 
symbols of this scheme are of the form 

aW(^)= ^f^' (27) 

lil<2n-l 

and since the scheme is interpolatory one has 

aW(z) + aW(-z) = 2, 2GC\{0}. (28) 

Due to the interpolatory definition of the subdivision it is clear that each function / G 
E (Aq..., A2n-i) is reproduced stepwise by the scheme, and Theorem [TO] implies that 

(-exp (-2-('^+i)A,)) = 0, . = 0,...,/x, - 1. (29) 

where /Xj is the multiplicity of Aj, i.e. the number of times the value A^ occurs in 
(Aq, A2n-i) • Hence the subdivision scheme based on interpolation in E {Xq, ...,X2n-i) 
is completely characterized by and (12^ . In the terminology of [13] this is the even- 
order, symmetric and minimal rank scheme reproducing -E (Aq, A2n-i) • Note that the 
Deslauriers-Dubuc scheme is a special case by taking Aq = ... = A2n-i = 0, reproducing 
the space of algebraic polynomials n2n-i- 
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Definition 17 For given real numbers Ao,...,A2n-i l^i a}^\k E No, be the symbols satis- 
fying ( f^) and ^''^y 'natural number m G Nq we define a new subdivision scheme 
S^" defined by the symbols 

aW'™ (z) = a['^+'"l (z) for k E Nq. 

In the rest of the paper we shall use the notations given in Definition [T71 

Remark 18 It is easy to see that the subdivision scheme S^^" is again an even-order, 
symmetric and minimal rank scheme reproducing E A2n-i/2'"). 

Many properties of the subdivision scheme S^" for exponential polynomials can be 
derived from their polynomial counterpart, the Deslauriers-Dubuc scheme. The key to 
these results depend on the following observation in [TU Theorem 2.7]. 

Proposition 19 The subdivision scheme S^" is exponentially asymptotically equivalent 
to the 2n-point Deslauriers-Dubuc subdivision scheme, i.e. there exists C > such that 

Yl \Pi - ^ ^^-'^ for all k E No. (30) 

lil<2n-l 

Deslauriers and Dubuc showed in [9] that their scheme is C'^-convergent implying 
the existence of a basic limit function which will be denoted in the following by <l>^2n 
Furthermore, one can find sufficient criteria for the C^-convergence in [9]. 

According to Theorem 2.10 in [M] we have the following result 

Theorem 20 Let £ E Nq. // the Deslauriers-Dubuc subdivision scheme is C^-convergent 
then the subdivision scheme S^" is C^-convergent as well. 

According to the last theorem the subdivision scheme S^" has a basic limit function 
which will be denoted in the following by . 

A function / : M"' — t- C is called Lipschitz function of order a (or Holder function of 
order a) if there exists a number L > such that for all x, y G 

\f{x)-f{y)\<L\x-yr. 

The set of all Lipschitz functions of order a is denoted by Lip{a). Functions in Lip [a) 
can be characterized via Fourier transform and we just recall Lemma 7.1 in P, p. 56]: 

Lemma 21 Let f : W ^ C be an integrable function whose Fourier transform is g {^) ■ 
We assume that l^f^'^g (0 integrable where £ G No and a E [0, 1] . If a = then f is i 
times continuously differentiable. If a ^ 0, then f^^^ is a Lipschitz function of order a. 
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Theorem 22 Let a G [0, 1) and £ G Nq and assume that the basic limit function $^2n qJ 
the 2n-point Deslauriers-Dubuc scheme satisfies for some e > and C > the inequality 



-l—a—e 



(31) 



for alluj G M. Then the basic limit function of the scheme S^" has its i—th derivative 
in Lip{a). 

Proof. We begin with some general remarks: let us define z^^^ = exp (— 2^*^'^+^)Aj) . 
Using (l29l) we can write 



aW(.)= n 



2n , \k] ' 



(32) 



Let D2n{z) = (^^)^" &z)2n(^) be the symbol of the Deslauriers-Dubuc scheme. It was 
shown in [131 Formula (2.32)] that fo^^' is exponentially asymptotically equivalent to &z)2„; 
so there exists a constant B > such that 



I^W (e-) -6^^^ (e^-)l <B-2 



-k 



(33) 



for all G M and for all A; G Nq. As an intermediate step we consider the non-stationary 
scheme S!^ defined by the symbols 



ctHz) 



z + l 



2n 



(34) 



According to the proof of Theorem 2.10 in [14j the scheme S!^ has a basic limit function 
denoted by Let be the basic limit function of 5*^^". By Proposition [7] 



1 



m+k—l]f iuj2 '° ^ 



(35) 



k=l 



Similarly, $^2n(^) = ]Tk=i ^o,{e"^'^ and we see that 



n 



n 

k=l 



Au]2 



(36) 
(37) 



It is well known ^ that the trigonometric polynomial fe^ij^ ( e ^"^^ 
the unit circle, hence the denominator in fl371) satisfies 



> 5 > 



does not vanish on 



(38) 
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for some S > 0. Using (155]) and fl33l) it is straightforward to prove that the infinite product 
in f l36p is uniformly bounded for w G M, and we obtain 



< Mr, 



Recalling that each factor induces convolution with a 5— spline of order (denoted 



as usually by Bq) the above inequality can be written as 

2n oo 



Comparing 
each factor 



and 



sm 



rn+k-l] ( ^iu2-'' 



2 ' k=l 

<M„C(|a;| + l)- 



-1—a—e 



(39) 
(40) 



we see that in order to prove the theorem we need to replace 
by the Fourier transform of the basic limit function generated by the 



subdivision scheme with symbols 



. This should be done with care since 



vanishes at infinite number of points on M. We employ another observation from [16 



that the basic limit function $o of the scheme with symbols 
-B— spline of order 0, 



z+z 



[k] 



is an exponential 







for X G [0, 1] 
otherwise 



Adding such factors, 



z+z 



[k] 



2n 



to the symbols {c''^'^} results in repeated convolutions 



of $^ with {Bq}'^,'!^. Since Bq(x) decays as A for \uj\ — y oo, each convolution adds one 



power to the decay power a in (1391) . Hence, after 2n convolutions we obtain 

^ 2 / k=l 

0((|a;| + l)- 



/ • \ 2n oo 



(41) 



sm 2 



A2„ 



It follows by Lemma ED that the (2?t, + £) — th derivative of is Lip{a). Now we are 
ready to return to by removing the factors 



from "^rnioj)- We need to show 



that each factor removed implies that the order of the derivative which is in Lip{a) is 
reduced by one. To show this we consider g = Bq * f , where / is a function of compact 
support and g' is in Lip{a). It follows that 

^'(^)=/(^)-/(^-l)- 
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Summing the above relations over all numbers s G {t + j}f=Q, for large enough N, we 
obtain 

N 
j=0 

implying that / is Lip(a). Since g{u) = ^^^/(w), we have just shown that the conse- 
quence of removing a factor is that the order of the derivative which is in Lip{a) 

2 

is reduced by one. Removing 2n such factors yields the desired result, namely, that the 
£—th. derivative of {^^"} is in Lip{a). m 

Remark 23 In JlBI, fSUf . f21^ motivated by the multivariate polyharmonic subdivision 
and wavelets on parallel hyperplanes, an explicit expression for the polynomial h^^^ {z) is 
found, for the case of frequencies given by Xj = ^, j = 0, 1, n — 1, for some real ^ > 0, 
and Xj = — ^, forj = n, 2n. The classical Deslaurier-Dubuc case corresponds to ^ = 0. 

4 Non- stationary multiresolutional analysis 

The concept of a multi-resolutional analysis, introduced by S. Mallat and Y. Meyer, is an 
effective tool to construct wavelets in a simple way from a given scaling function see 
e.g. [I], [H]- Non-stationary multiresolution analysis was introduced in [8] by C. de Boor, 
R. DeVore and A. Ron. For convenience of the reader we recall here the definition for the 
univariate case: 

Definition 24 A (non-stationary) multiresolution analysis consists of a sequence of closed 
subspaces Vm,m G Z, in (M) satisfying 
(i) Vm C Vm+1 for all m G Z, 

(a) the intersection flmg^Kn is the trivial subspace {0} , 
(Hi) the union UmezKn is dense in (R), 

(iv) for each m G Z there exists a function ip^ G Vm such that the family of functions 
{ipm (2™t — A;) : A; G Z} form a Riesz basis of Vm- 

The function (pm in condition (iv) is called a scaling function for Vm- The requirement 
(iv) means that for each f & Vm there exists a unique sequence (cfc)^^^ in (Z) (i.e. that 
\ckf < oo) such that 

oo 

f{t)= Yl Cfc<^™ (2™t - A;) 

fe=— oo 

with convergence in (M) and 

oo 
fc=— oo 
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J2 Ck^m (2"t - k) 



< Br, 



E 



for all (cA;)^g2 in P (Z) with < < -Bm < oo constants independent of / G Vm- 

Using (i) one can define the wavelet space Wm as the unique subspace such that 
Vm © Wm = Vm+i for m G Z and Wm is orthogonal to Vm- It easy to see that this implies 
that Wk and Wm are orthogonal subspaces for k ^ m. Conditions (ii) and (iii) imply that 

(R) = ®mezWm. 

We refer to [8] for an extensive discussion on the construction of so-called pre-wavelets 
for a nonstationary multiresolutional analysis, see also Important examples of non- 
stationary multiresolutions occur in the context of cardinal exponential-splines wavelets 
which generalizes the work of CK. Chui and J.Z. Wang about cardinal spline wavelets in 
[4J, see [2], [3]. The interested reader may consult [29], [32] for the theory of exponential 
splines and [8], [19], [23], [25], [28] for the construction of wavelets in this context. In 
passing we mention that the results in the last cited papers have been rediscovered in [35] . 

Definition 25 A multiresolutional analysis is called stationary if in condition (iv) the 
scaling function (pm &V is the same for all m G Z. 

In this paper we shall be concerned only with orthonormal non-stationary MRA: 

Definition 26 A multiresolutional analysis is called orthonormal if in condition (iv) the 
functions t i — > 2™/^™ (2^"^ - k) for k e Z are an orthonormal basis of Vm- 

Let us recall that an orthonormal wavelet is a. function in (M) such that the 
system ^pm,k (x) = 2™/^^/' (2™a; — k) with m, /c G Z is an orthonormal basis of (R) . In 
the context of nonstationary wavelet analysis one wants to find a sequence of functions 
1pm e (R) , m G Z, such that 

lljm,k (X) = T'l^lpm (2™X - k) 

with m. A; G Z is an orthonormal basis of (R) . 
5 Daubechies type wavelets 

Daubechies wavelets tp are orthonormal wavelets with compact support and certain degree 
of smoothness. The construction of Daubechies wavelets is often presented in the following 
way (see e.g. [1], [T7], [27]) : using the concept of an orthonormal MRA it suffices to 
construct a suitable scaling function y?. Elementary considerations show that the Fourier 
transform ^ of the scaling function Lp should be of the form 

oo 
k=l 

where m {u) is a trigonometric polynomial with real coefficients and m (0) = 1 satisfying 
the equation 

\miuj)f + \m{uj + n)f = 1- (42) 
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This leads to the question which non-negative trigonometric polynomials q (u) satisfy an 
equation of the type 

q{u) + q{u + n) = l and q (0) = 1. (43) 

There are many explicit solutions of (H3|) . For example, if n is a natural number then the 
trigonometric polynomial 

qn (to) = 1 — c„ y (sint)^"^^ dt 



with c„ := f sin^"~^ satisfies equation (1431) . By the Fejer-Riesz lemma one can find a 



(non-unique) trigonometric polynomial m (u) such that 

q:,{uj) = \m{u)\\ (44) 

We call a Laurent polynomial m {u) with real coefficients and m (0) = 1 satisfying fH^ a 
Dauhechies filter of order n. The Daubechies scaling function ip"^ for the Daubechies filter 
m {oj) is then defined by 

oo 

(^(w) = JJm(2-^u;) . 

k=l 

This procedure is elegant but the construction of the trigonometric polynomial g„ in the 
above approach seems to be rather miraculous. Let us emphasize that I. Daubechies has 
shown much more (see e.g. [6[ p. 210] or [37]): the regularity of the wavelet and the 
scaling function imply that the symbol m {u) must contain a factor (1 + e*"^)" /2". Hence 
g„ (w) is of the form 

(1 + e*^)^" 

^" = ^ ^2n-l (W) 

where -F2„_i (w) is a suitable trigonometric polynomial with real coefficients which can be 
determined by Bezout's theorem from (H3|) . Indeed, it follows from these considerations 
that 

where is the symbol of the Deslauries-Dubuc subdivision scheme, a fact which is 
already mentioned by Daubechies in her book [6, Section 6.5] giving credit to this obser- 
vation to M.J. Shensa in [33], see [6l p. 210]. Hence the Deslauriers-Dubuc scheme leads 
in a very natural and direct way to the construction of the Daubechies scaling function 
and therefore, by MRA-methods, to Daubechies wavelets. 

In this section we want to use this concept in order to define Daubechies type wavelets 
for exponential polynomials. In this setting we have some additional freedom which 
is interesting for applications: we may choose real numbers Ao,...,A„_i and we shall 
construct Daubechies type wavelets reconstructing the space E (Aq, A„_i) . In the case of 
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Daubechies wavelets this corresponds to the fact that the Daubechies wavelet reproduces 
polynomials of degree < n — 1. . 



We shall write shortly Aq = (Aq, A„_i) and define Xn+j '■= —Xj for j = 
1. We consider now the subdivision scheme based on interpolation in E (Aq, 
According to Definition [T7I the subdivision scheme Sq'^" has the symbols 



2n- 



n 



1 , [k] 



with zf^ := exp {-2-^''+^^ Xj) 



0, ...,n - 

■■, X2n-l) ■ 



(45) 



for k E No and j = 0, ...,2n — 1. Crucial is the result by Micchelli in [30] (Proposition 
5.1), proving that the symbols a^'^l satisfy 

a[^'l(z) > , for \z\ = l , (46) 

with equality possible only if z = —1. By the Fejer-Riesz lemma we can find for each level 
k a "square root" Laurent polynomial M^''^{z) with real coefficients, satisfying 



Note that this implies that 



-|MW(e^'^)P and M^''^ (1) > 0. 



(47) 



,[k] 



-MW(z)M[^1(-) 



2 -z- 

for all complex 2; 7^ 0. Again, there are many Laurent polynomials M'^^^ [z) which satisfy 
pri) and all possible choices can be described through suitable subsets of the zero-set 
of a}'^\z). First we choose the roots z = — exp (— Aj/2'^+-'^) for j = 0, ...,n — 1, in order 
to obtain stepwise reproduction of the space E {Xq, Xn-i) , see Proposition U7\ below. 
Further, we have to choose another n — 1 roots of the factor fol'^l in fH^ . Since IS 
symmetric, its 2n — 2 roots come in inverse pairs, say Zi and 2;"^, and as well complex 
conjugates and 11'^ if Zi is not real, for i in an index set In-i- We choose either the set 
{zi,zl} or the set \^z~^ ,zi~^^ for each i G /„_i, leading to a Laurent polynomial with real 

coefficients which still has to be normalized so that M^'^^ (1) = A/2aW(l) > 0. We shall 
call a sequence of filters M''^' (z) , k E Nq, chosen in this way a non-stationary Daubechies 
type subdivision scheme of order n. 



Since Ml^l (1) is positive it follows that 1 < iiVfl'^l 
infer that 



1) + 1 and aW(l) = iMW(l) 



we 



-Mt^'] (1) - 1 
2 ^ ' 



< 



-AfW (1) - 1 



2 



-aW (1) 
2 ^ ' 



Since a^^l (z) is exponentially asymptotically equivalent to the Deslauriers-Dubuc scheme 
and hD2n (1) = 1 we infer that there exists C > such that for all k E Nq 



-MW (1) - 1 
2 ^ ^ 



At first we notice the following result: 



<C-2-'' 



(4J 
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Proposition 27 Let Ao,...,A„_i be real numbers. Then there exists G Nq such that 
the Daubechies type subdivision scheme reproduces stepwise functions in (Aq, A„_i) 
for all levels k > ko. 



Proof. Let z^!'^ = exp (— Aj/2'^"'"-'^) . By construction M^^^ {z) has a zero at —Zj''' of 



multiphcity fij, the number of times Xj occurs in (Aq, A„_i) , hence 
dx^ 



(^-^r) = for s = 0, - 1 and j = 1, ...,n- 1. (49) 

By (H71) and the fact that jat*^]} reproduces stepwise functions in E {Xq, X2n-i) we 
conclude that 

-|MW(zf])|2 = aW(4']) = 2. 
2 

Since zj'^^ is real and iVft'^l (2;) has real coefficients it follows that Mt^l(2;j^') is real so 
j\^[fc]j^^W-j = 2 or —2. Since zj'^^ converges to 1 for k ^ 00 and M[^l(2;j^') converges to 
M^2n(i) > there exists ko G No such that MW(^].^') > for all k > ko and j = 1, ...n-1. 
Hence M['^](z'-^^) = 2 for all k > ko- From fHTl) we infer that for real x 



For s = 1 this means that 



-^aW (x) = y r -^mw 



r=0 



r / fix'' 



[k] 



+ 



dx 



I z 



[k] 



M'*^! has real coefficients and is real we conclude that -^M^^^ ( z 



dx 



[k] 



we obtain that -^M^^^ 



[k] 



for s 



1. 



[k\ 



. Since 



0. Inductively 



Proposition 28 The product ImI^^-^I | e 2^ 



converges. 



Proof. By construction M^^^ (e**^) has real coefficients and M^^^ (1) > 0. Since 

< 2, and therefore \Mk < 2. Moreover, it follows 



we infer from (EHD that \a^^^ {e'"^) 



from (gHD that 



E 

fc=i 



-MW (1) 



< 



c5;]2-^ 



A:=l 



Proposition [6] finishes the proof. ■ 

In order to obtain asymptotic equivalence for the non- stationary Daubechies subdivi- 
sion scheme we have to choose the filter M^^^ iz) with more care: 



Theorem 29 Let M^^" be the Daubechies filter of order n as defined above and let 
Ao,...,An-i be real numbers. Then there exists a non-stationary Daubechies type sub- 
divison scheme {M^''^(z)} which is exponentially asymptotically equivalent to M^^" and 
reproduces e^^^ stepwise for all k > ko and for j = 0, n — 1. 
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Proof. Note that M^^n has n zeros at —1 and n — 1 other zeros, say ai, 
which are of course zeros of the factor Bd^^ (z) of the Deslauriers-Dubuc symbol (z) = 
^i^-j2n^^^^ ^ Using (!24|) and Proposition [T6] we see that ai, are pairwise dif- 

ferent and simple zeros of (z) . Recall that at'^l (z) is exponentially asymptotically 
equivalent to the symbols (z) . Then z'^"'~^a^''^ (z) are polynomials and z'^'^~^a^''^ (z) is 
obviously exponentially asymptotically equivalent to z^"~^Z)2n {z) . By Theorem [T5] there 



^2n 

exists a constant C > and a zero af^ of at'^' (2;) such that 



[k] 

a J — aj 



< 02'" for all 



k gNq and j = 1, n — 1. Take ko G No large enough so that: 

(i) for each k > k^ the balls \z ~ OLj\ < C2~^ have empty intersection with the unit 
circle, 

(ii) they are pairwise disjoint for j = 1, ...n — 1, and 

(iii) they have empty intersection with the x-axis if aj is a non-real zero. 
Then for each k > ko there is for given a,- exactly one zero with 



[k] 

— aj 



< C2-*^ for J = l,...,n- 1, (50) 



leading to a unique choice for M''^' (z) for k > ko- Further the leading coefficient c'^^^ of 
the polynomial z'^^'^M^''^ (z) is determined by the equation 



71—1 71—1 



mw (1) = cw n (1 + ?0 ■ n (1 - ) • (51) 

i=o i=i 

Let c be the leading coefficient of ^^"^^£'2^ (^)- By (HH]) and ({50]) and (EI]) it is easy to 
see that there exists D > such that 

|cW -c| < L)2-'= 

for all k E Nq. By Theorem [T3] the subdivision scheme defined by the symbols z^"~^M['^1 (z), 
/c e No, is exponentially asymptotically equivalent to the scheme defined by the symbol 
^2n-ij|^_D2n_ This implies that M'^^ (z), k E Nq, is exponentially asymptotically equivalent 

to M^2n_ g 

Remark 30 Assume that M^^" is the Daubechies filter such that all zeros 7^—1 have 
absolute value bigger than 1. Then one can define (z) in the last theorem by the 
condition that all its non-trivial zeros have absolute value bigger than 1. 



It thus follows, from the theory of asymptotically equivalent schemes in [TB], that 
the scheme with symbols {M^'^'^™'^ (z) , k E No} defines continuous basic limit functions 
{v'm (■)}• Proposition [7| shows that 



e 2* 



00 

j^[m+k-l] 

^5M=n — ^ — -■ (52) 



k=i 
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In particular, we have 



2 



{uj) (53) 



where is the basic hmit function of 5*4^". 

The following observation is straightforward and the proof is included for convenience 
of the reader: 

Proposition 31 The functions (2™ ■ -k),k E Z, are orthonormal. 



Proof. Define Dk^i := {T^t - k) Lpt° (2™t - l)dt. A simple transformation 



of variables and Plancherel's formula gives 



1 



Dk,i = / vt' {y - k) {y - l)dy = — / {u) {u)e-'^'-'^-dt. 

J -oo ■^'^ J -oo 

Using (l53ll we obtain that 



Definition 32 For each m E Nq we define the linear spaces Vm by 

' V oo }> . (54) 



Vm:=\fe L'iR) I fit) = Y.c,^^^{2-t-3) , l^^'l 



We remark that we could also define Vm for integers m G Z since the symbols 
and the scaling function and ip!^ could be defined for all m G Z. 

Proposition 33 The spaces Vm are nested, i.e. that Vm C Kn+i for all m G Nq. 

Proof. Using the product representation (1521) we obtain 

(c) = n = ^^^A^ (I) • (55) 

k=l 

Let us write M^"^^{z) = X]j=-n+i f^^P^^'' where fJj^^ are real numbers. Using elementary 
techniques in Fourier analysis it is easy to see that the equation fl55|l is equivalent to the 
refinement equation 

n 

¥>t'it)= E /^rV™V(2t + j)- (56) 

j=-n+l 
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Replacing t by in ( ISBj) we obtain that (p^{2"^-) G Ki+i- Similarly it follows that 
ifi^,^i2^t - j) E Kn+i for each j G Z. ■ 

We note that the orthonormality of ip^' {t — I) implies that 

/OO 1 
Vt' it) {t - l)dt =-J2 f^'^f^^i- (57) 

Indeed, using in the integral in and then again the orthonormality relations one 
obtains 



/OO 
^mV(2t + J>:^%(2t-2/ + fc)cit 

1 fOO 

= ^ E z^^/^lr^ / - J - 2/ + fc)rfy. 



Proposition 34 The union of the spaces is dense in L2 (M) . 

Proof. Define fm (x) = ift^ (2™a;) . Then £ (w) = ipt (2-'"a;) . By Theorem 4.3 in [8] 
it suffices to show that the set Q, defined as the union of the supports of fm over m G No, 
is equal to M minus a set of Lebesgue measure 0. Since $m"(w) = (w)P we have only 
to investigate the zeros of (u). By (!52|) this is an infinite product of non- negative 
terms a^^'^^^ ^e*"^^ j which may be zero only if e*'^^ = — 1, see [30]. Hence, the zeros 

of $m " are of the form u = 2'^ {tt + 2n7c) for some n E 1^. Hence the support of $m " is 
equal to the real line. This ends the proof. ■ 

Theorem 35 Let M^'^"{z) be a Daubechies filter of order n and assume that Ml^l (z) is 
as in Theorem Let a G [0, 1) and i G No, and assume that the scaling function ip^^" 
of Daubechies defined by the symbol M^^"{z), satisfies for some e > and C > the 
inequality 

<c(|w| + i)-^-i-"-" 

for alluj G M. Then the scaling function y?^" associated to the subdivision scheme M^^^™-^ {z) 
A; G No, have i~th derivative in Lip{a). 



as 



Proof. Let us define Zj = exp (^_2-{k+i) x-) . Then the symbol {z) can be written 
Similarly we can write for the Daubechies filter 



n 
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Since Ml*^! (z) is exponentially asymptotically equivalent to the Dauchechies filter M^^"- 
and B^''^{z) has only simple zeros it follows from Theorems \T5\ and [T^ that -B'^^ is expo- 
nentially asymptotically equivalent to B^'^'^{z), so there exists a constant C > such 
that 

Now one can proceed as in Theorem [221 ■ 

Next we turn to the construction of Daubechies type wavelets. The procedure will 
follow the classical pattern in MRA. For convenience of the reader we shall briefly sketch 
the construction: Recall that M^'^^{z) = Y^=-n+i I^^T^ write the refinement 
equation in (1561) in the form used in MRA, namely 

The Daubechies type wavelets tp^ are now defined in the classical way, namely by 

^^"w = E^fV:^'Vi(2t-j), (58) 

where the coefficients are related to those in f l5B]) by 

= (-1)^+VL% . (59) 

Then il)^ has compact support since i^^j^i has compact support. It is a routine exercise 
to see that the system of functions \^2^^'^-ip^ {T^t — r) : r G Z} is orthonormal: define 

oo 



Dr,s := 2"^ / (2™t - r) (2™t - s)dt. 



The transformation y = 2"^t — r gives Dr,s = j"^^ i^m iv) {y + r — s)dy and therefore 

/oo 
^t^A'^y-k)^t^,{2y + 2T-2s-l)dy 

1 roo 

= o E 4"^'^!"^ / {x + 2r-2s + k- l)dt 



1 X ^ [m] [m] 1 \ ^ [m] [in] i- 

k&Z k€Z 

where we have used equation ( 157]) . 

As explained in Section 4 the wavelet space Wm is defined as the orthogonal comple- 
ment Wm of Vm in Vm+1- Ncxt one shows that the wavelet space Wm is equal to 

Wm:=lfe L\R) I fit) = Y,c,2^l''^t{'^^t - 3) , |c,f < oo I . (60) 
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We show at first that Wm is orthogonal to C Kn+i implying that Wm C Wm- So we 
look at 



oo 



Cr,s := 2"^ / ift' (2™t - r) iji'{2^t - s)dt. 



Again Cr,s = JZo iv) + ^ " ^)^^ ^nd (|56D and dSHD yield 

/oo 
^:l%(22/-J>:^%(2y + 2r-2.-/)ciy 

1 /-oo 

= ^ E ^'^-M"'^ / ^™Vi(^)v':^°+i(^ + J + 2r - 2. - 

A simple well known argument shows that the last sum is always zero, see e.g. [271 P- 123]. 
The proof that Wm © is equal to Vm+i (implying that Wm = Wm) follows standard 
arguments in MRA and is omitted. 

Corollary 36 The smoothness of the Dauhechies type wavelets ipm least as that of 

the classical Dauhechies wavelet ip 

Proof. This follows immediately from formula ( l58ll and Proposition |35l ■ 
Finally we mention that the concept of reproduction is defined in MRA in the following 
way, see e.g. [36] : 

Definition 37 A non-stationary multiresolutional analysis {ym)m&z '"'^^^ compactly sup- 
ported scaling functions ipm reproduces a function / : M — )■ C if for each m G Z there exist 
complex coefficients Cm such that 

f{x)=Y,Cn.Vrn{'2"'x-l). (61) 

Note that the sum in (16T1) is well defined since is compactly supported. It is proved 
in [36] that fH9|) implies that the MRA {Vm)mez Dauhechies subdivision scheme 

reproduces the space E (Aq, A„_i) and the wavelets ipm have vanishing moments in the 
sense that 

j (t) e^^^dt = for all j = 0, ...n - 1. 

As we mentioned in the introduction an attempt to construct (non-stationary) Dauhechies 
type wavelets based on the subdivision scheme for exponential polynomials in analogy to 
the case of Deslauriers and Dubuc, was carried out in [36]. The authors showed that 
many results and techniques of classical multiresolutional analysis, briefly MRA, carry 
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over to non- stationary MRA (except that all filters and scaling functions are dependent 
on the scale of the multiresolution) and they introduced non-stationary Daubechies-type 
wavelets reproducing a family of exponentials {e'^J*}_^._^ under the assumption that the 
filters have sufficiently long support. This strong assumption has been made in order to 
guarantee that the symbols of the nonstationary scheme are non-negative. However, as 
we mentioned in the introduction, a main result proven in [30] shows that the symbols 
are non-negative under the simple assumption that the exponents for j = 1, ...,n are 
real. As a consequence, the question about the order of regularity of the new Daubechies 
type wavelets, is not addressed in [36j- 
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